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ABSTRACT 

Magnetic Resonance Imaging (MRI) techniques are used to scan images of brain tumors. The extracted MRI image is sensitive to 
noise which limits the visibility of certain characteristics of the image. This unwanted noise can be disintegrated and curtailed from 
the original Image by using the non linear digital median filter. The intensification is achieved by using the combination of discrete 
wavelet transform (DWT) and stationary wavelet transform (SWT) augmented with Bi-cubic interpolation algorithm. The quality of 
the resolute image is analyzed with the possible static parameters. Segmentation and clustering approaches are applied on the 
resultant sharpened high resolution image to estimate the tumor area. In this paper the results obtained with proposed Optimized 
Hybrid Clustering (OHC) algorithm is compared with the existing K-Means, Pillar K-Means and Fuzzy C-means (FCM) approaches. 
The tumor sensitivity and accuracy of the test parameters with the OHC approach are estimated. In-order to deliver the drug to kill 
the aberrant tissue the number of tumor cells are estimated against its radius. The Proposed method exhibits reliable values when 
compared with existing approaches. 
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1. INTRODUCTION 

The cells in a human brain will develop in a controlled manner. The old and flawed cells will be recovered with new cells. The cell 
division grows in undisciplined manner and form as an aberrant tissue labeled as Tumor. Tumors are classified as primary and 
secondary tumors. Treatment of the tumors depends on its type, area and location. Magnetic Resonance Imaging (MRI) is adopted 
as reliable and superior approach to trace the brain tumors. MRI Scan is a non invasive procedure that uses strong magnetic field of 
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21.1 Tesla augmenting with radio frequency waves as pulses to attain the images of the brain. Three gradient magnets used in MRI 
labeled as X, Y, Z. Each gradient field aligned at different plane of the organ. The gradient magnetic field focuses on the distinct area 
of the slice and work in association with radio frequency pulses to produce the image. The sequence of X, Y, Z gradients used to 
generate image of the slices in any orientation. This is one dominant component forced to adopt MRI as a diagnostic tool.MRI 
image consists of Grey and white pixels. The grey pixels represent both signal intensity and noise. By increasing the electrical 
currents in the wires of the gradient magnetic field the main magnetic field opposes it and generates huge noise. And irregular 
measurement of electronics also introduces significant noise. This noise minimizes the visibility of certain characteristics of the 
image. Noisy image may introduce inadequate allotment of Gray level pixels 

as white level pixels. The visibility of the low contrast MRI images can be improved by removing the noise with a median filter. A 
Non linear digital median filter used to remove the noise in MRI images and safeguard the edges of the image. Each pixel in the 
image relates with nearby neighbours to determine the representative of its surroundings. The median of neibouring pixel value is 
determined to modify the pixel values. The median value will not be affected with aberrant single pixel in a neighbourhood. So the 
median is more potent average than mean. Using this approach to save sharp edges in MRI brain tumor images, the filter does not 
generate improbable pixel values. The substantial noise is withdrawn from the MRI brain image with median filter [1][2]. Image 
enhancement in MRI images is used to resolve the problem of low resolution and to increase the contrast and sharpness of the 
image. The de noised image is further enhanced using combination of discrete wavelet transform (DWT) and stationary wavelet 
transform (SWT) with applied bicubic interpolation approach to generate sharp high resolution images. The high frequency 
components of the MRI image are perpetuated using stationary wavelet transform method. The loss of report arises by virtue of 
down sampling in discrete wavelet transform (DWT) begin in respective to sub bands. SWT is used to curtail the information loss. 
The discrete stationary wavelet transform disintegrate the MRI image into low and high frequency sub bands. [3][4][5][6] Bicubic 
interpolation is then applied to sub band images. Bicubic interpolation is achieved using nearest neighbour image resampling. The 
re-sample images obtained are stable. The proposed SWT with Bicubic interpolation method generates sharper high resolution MRI 
image [7]. The enhanced high resolution MRI image quality is analyzed with static parameters. The quantitative analysis concludes 
that the combination of DWT+SWT with applied bi-cubic interpolation results reliable resolute image. Which was further processed 
using segmentation and clustering algorithms. The Optimized hybrid clustering (OHC) algorithm enforced to cluster the tumor 
associated pixels. These pixel values are used as optimum values to estimate the area of the tumor from MRI images. 


2. TECHNICAL REVIEW 
2.1 De-noising 
Transform domain and spatial filtering approaches are adopted for MRI images. The compelling noise in MRI images is curtailed by 
adopting spatial filters. Mean and wiener filters are linear and median filter is nonlinear. The salt and pepper noise in MRI images is 
curtailed with non linear digital median filter. While removing the noise the filter safeguards the edges. An array of MRI Pixel points 
are glanced to elect the representative of its surroundings. The median of neighboring pixel value is to be determined to refine each 
pixel value of the MRI image. The median filter is alleged to be edge secure. The estimated median value will substitute the pixel 
examined with the center pixel value. If the considered neighbourhood has an even number of pixel points, the average of the two 
centre pixels is used. 

TG. y) =1) 
I, <I, = I, _— a 


Where, T is the threshold value and {a * m/2)is the centre pixel intensity value. As edges contain important information for 
segmentation, median filter is used for preserving significant details in MRI image. Using wiener filter compelling noise in an image 
is to be diminished by correlate with desired noise free image. Mean wiener and median filters are used for comparative study of 
MRI brain tumor images. Wiener filters update the pixel value by averaging of its neighborhoods. 


The Maximum signal power to the maximum signal noise is indicated as peak signal-to-noise ratio (PSNR). The PSNR is determined 
using the following equation 
PSNR ab =10 Log 10 255*/ MSE 
Where, MSE = aad dy aallGr, y) —I'@ y)F* 
Where , M = Number of rows in an image 
N= Number of columns in an image 
MSE is mean squared error, | is an input image and I is output image [1][2]. 


2.2. Image intensification Ne) 
Discrete wavelet transform and Stationary wavelet transform are used to resolute the Low resolution MRI input image. The ‘© 
information loss on its edges i.e high frequency components may arise by applying the interpolation for image intensification. 
Safeguard the edges is imperative. DWT has been hired to safeguard the high frequency components of the image. In-order to 
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obtain the redundancy and shift invariance the discrete wavelet transform coefficients are essentially interpolated. Hence DWT is 
used to disintegrate the MRI image into sub band images. The high frequency components of the image labeled as sub bands LH, 
HL, HH. 

For interpolating data points on a two dimensional regular grid, Bi cubic interpolation is used in order to obtain smooth surface. 
A growth factor of ‘2’ for Bi cubic interpolation is enforced to high frequency sub band images. Stationary wavelet transform (SWT) 
is adapted to negotiate the information loss in the Sub bands due to down sampling. The MRI brain tumor image disintegrates into 
different sub band images by using stationary wavelet transform, low-low (LL)), low-high (LH), high-low (HL), and high-high (HH). 
The SWT High frequency sub bands and interpolated high frequency sub bands are combined each other. Further intensification is 
achieved by interpolating the estimated high frequency sub bands. The authentic high resolution input image is having admirable 
information than low frequency sub bands. Hence the interpolation is again applied on the input MRI image Low frequency sub 
bands [3][4][5][6]. 


2.3. Segmentation and Clustering 

Image segmentation and clustering practice are imported to estimate the area of the tumor. Image segmentation is classified into 
Pixel based, regional and edge based methods. In this paper the brain tumor images are sub divided into multiple segments as sets 
of pixels using pixel based segmentation. The MRI image incorporates white and grey color pixel elements. White color pixel data 
points are related to tumor cells and the Gray color pixel data points relate to normal cells. Collection of data point of the pixels that 
belongs to the same color will be quantified using Euclidian distance method. The clusters may contain large number of pixels. The 
pixels may be either close or far from the cluster centre. 

Define the cluster centers and each point is to be set aside to the closest cluster centre. Then every centre is the mean of the pixel 
points allocated to that cluster. In order to estimate the area of the tumor “Optimized hybrid clustering algorithm” is applied on the 
Input MRI images. 


a) Feature Extraction 

A Binary Mask is applied on the MRI input image. Grey and white pixels become brightened. A threshold value is defined to 
compare with the coefficients of each pixel. If the value lies within the threshold value a ‘Zero’ is assigned to that coefficient else a 
‘One’ is assigned. The Fuzzy C Means output is the extracted tumor cluster from the MRI image. The magnitude of the coefficients 
from the extracted tumor cluster is above the threshold value. The threshold value is determined using the MATLAB Command 
(Global Image Threshold Function) the threshold value is represented as ‘T’. 


T= Max [ f(n,) ] + Min [f(n)] / Avg of f (nk ) 


The MRI input image is represented as ‘f’. Grey level pixels are represented with ‘k’. Each pixel in an image ‘f’ is compared with the 
threshold value T. The threshold value lies in the grey scale range of ‘kK’. In-order to define a value to the specific pixel from the 
output image ‘g’ a binary decision is applied.[7][8][9]. 
g (n) = '0' if f (n) > T 
='l' if f(n) <T 


b) Approximate Reasoning 
Linearization method is adopted to estimate the area of the tumor. The MRI image comprises of two pixel values either grey or 
white (0 or 1). The size of the image is considered to be 256 % 256. The aggregate of white and black pixels represents the Binary 


image. [7][8][9]. 


255 255 


1=> > F@+Fap 


==o0 y=0 

Pixels= Width * Height = 256*256 

f(0) = black pixel ‘0’ 

f(1)= White pixel ‘1’ 

Where, 

Totel White Pixels (4%) = fcay "= number of white pixels 
i ee 1 Pixel value is equated to 0.264 mm 

The equation to estimate the tumor area is 


Tumor Size (S) = (|w, i 0.264) 


3. PROPOSED METHODOLOGY 

Optimized Hybrid Clustering (OHC) algorithm is a combination of pillar k-means and Fuzzy C means approaches. The combination 
of these algorithms favorably reduces the difficulties occurred in previous segmenting algorithms. The membership values to the 
pixels are considered -1, 0, and 1 using cluster centroids cause to minimize grey level similarity limitation with FCM alone. The 
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floating points obtained with K-Means due to mean calculation for new centroid updating in some iterations has been favorably 
curtailed by considering the maximum data set point value. The OHC Algorithm substantiates the implication of MRI brain tumor 
image segmentation and clustering to estimate the tumor area. This approach endorsed optimized results. 


3.1 Optimized Hybrid Clustering Algorithm 

Step 1. Normalize the Input image into 256*256 sizes. 

Step 2. Convert the image into double precision values. 

Step 3. Reshape the image into data sets of column vector for segmentation. 

Step 4. Determine the number of clusters i.e. centroids. 

Step 5. Calculate the distance between input data set points and set points of the centroids. 

Step 6. Find the data set points which are neighboring to the centroid. 

Step 7. Select the centroid with minimum distance then move the data set points to the closest relevant centroids 
Step 8. Re compute the centroids by selecting the maximum pixel value from the set of relevant centroid data points. 
Step 9. Repeat the process until the new centroids and the previous centroids are symmetrical. 

Step 10. Allocate the membership values -1, 0 and 1 to each pixel of the cluster centroids 

Step 11. Estimate the global threshold value T. 

Step 12. Determine the data set points with similar membership value and then reform the new clusters. 

Step 13.Apply the binary decision to obtain the resultant cluster image g (n) 

Step 14. Estimate the square of number of white pixels (Wp ) 

Step 15. Compute the area of the tumor “S” using linearization method 


4. RESULTS & ANALYSIS 

An innovative procedures based on Median Filter, wiener filter and mean filters Algorithms are adopted for effectively de noise the 
MRI brain tumor images and for preserving the edges. While preserving high quality of restored image the noise abolition across an 
ample range of noise quantity about 30 to 98% is considered. This is mainly due to median filter used for the MRI images. 
Considering significant number of MRI brain tumor images, Wiener filter deteriorated to preserve the fine point pixel values and 
edges of the MRI image. The PSNR estimated values for the above three filters are tabulated. Table 1 and Fig.6 shows the median 
filter produced favorable PSNR Results when compared with other filters. The median filter applied to different MRI images caters 
satisfying results. The DWT and SWT used augmenting with Bi-cubic interpolation is enforced to enhance the de-noised tumor 
image. Testing with interpolation factor 2 applied on DWT and SWT to resolute the image influenced on the high frequency 
components. Considering fig 1(c), fig 2(c) introducing the SWT information loss due to down sampling in DWT is considerably 
reduced. DWT has been selected to preserve the high frequency components of the MRI image. The image size of 256*256 is used 
for enhancement with an interpolation factor of 2.The Qualitative parameters of the image are tabulated. Table 3 and fig 5 
summarizes the quality of the resolute images. The laplacian mean square error Results less than 1.2 for the tested samples. The 
normalized absolute error results ‘0.5’. Less than 1 considered as a best value. The maximum difference and average difference 
values are approximately closes to zero. The peak signal to noise ratio is also favorably good for all the tested samples. Low level 
structural information i.e within 3.0 when compared with input image and resolute image. The value of SSIM represents ‘0’ labeled 
as zero correlation and ‘1’ labeled as exact correlation. The table 4 concludes the SSIM values for all the tested samples are 0.99. le. 
99 percent of the exact correlation is achieved. Table 2 summarizes the estimated Area of the tumor for sample 1, sample 2. Manual 
segmentation and Fuzzy C-Means procedures produce unfavorable results. Pillar K-means produced the tumor area for sample 2 is 
14.7013mm2. But, the proposed “optimized hybrid clustering approach” proved itself as reliable result for sample 2 with 10.9679 
mm2. Table 2 and fig 5 concludes that OHC approach produces the optimum tumor area for all the tested samples. While 
comparing the fig1(d), 2(d), The OHC approach produced the output exactly matched with the actual tumor shape. The total white 
pixels are precisely detected causes to estimate the true area of the tumor. This approach achieved the precision level at about 98% 
closest to the true tumor size. 


Tested Results 
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Table 1 PSNR values for various filters 
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Fig.2 (b) De noised images 
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Figure 3 PSNR values for various filters 


Table 2 Area estimation for various algorithms 


: Area of the tumor( mm’2 ) 
.N Al h 
eNe cei Sample Sample2 
1 | Manual 7.7914 23.2185 
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Table 3 Quality analysis for various sample images 
Quality parameter Sample1 Sample2 
PSNR(db) 65.00784 69.202903 
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SC 3.644875 3.956706 
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Figure 5 Quality parameters for resolute images 


5. CONCLUSION 

The experiments were conducted on primary brain tumor MRI images. The median filter method is produced favorable PSNR value 
and justified as an optimum denoising method for MRI images. The combination of DWT-SWT in conjunction with bicubic 
interpolation approach results intensified quality image analyzed with static quality parameters. Based on the statistical parameters 
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the resolute image is evincing acceptable quality image. The proposed “OHC” approach proved as a refined and respectable method 
to estimate the absolute area of the tumor. This Proposed OHC algorithm can be applied to any MRI image for estimating the tumor 
area and volume. The number of tumor cells are estimated w.r.t its radius. Specifically the extreme dependence on MRI to detect the 
tumor cells, with the advances in tumor diagnosis, can increase the survival of tumor patients. Reckoning the tumor area, volume, 
and number of aberrant tissue cells will play a vital role to treat the tumors. Subsequently the location of the tumor may be 
determined in addition with the size i.e. area of the tumor. Precisely the Qualitative and quantitative analysis of the tumor will play a 
vital role for applying the radiation to annihilate the tumor cells. 
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